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Abstract 

The ZEUS inclusive differential cross-section data from HERA, for charged and 
neutral current processes taken with and e“ beams, together with differential 
cross-section data on inclusive jet production in e+p scattering and dijet produc¬ 
tion in 7 p scattering, have been used in a new NLO QCD analysis to extract the 
parton distribution functions of the proton. The input of jet-production data 
constrains the gluon and allows an accurate extraction of as{Mz) at NLO; 

as{Mz) = 0.1183 ± 0.0028(exp.) ± 0.0008(model). 

An additional uncertainty from the choice of scales is estimated as ±0.005. This 
is the first extraction of as{Mz) from HERA data alone. 
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1 Introduction 


Since the advent of HERA, considerable progress has been made in the determination of 
the parton distribntion fnnctions (PDFs) of the proton. Precise knowledge of the PDFs, 
and of the strong conpling constant, as{Mz), is crncial for an nnderstanding of proton 
structure. Moreover, it is required for any calculation of cross sections at hadron colliders 
both for Standard Model physics and for the discovery of physics beyond the Standard 
Model. 

The PDFs are usually determined in global hts [1-3] made within the conventional 
DGLAP formalism [4-7] at next-to-leading order (NLO). Such hts use data from many 
different experiments, with the inclusive cross-section data from deep inelastic scatter¬ 
ing (DIS) experiments providing the major source of information. The wide kinematic 
range covered by the HERA DIS data [8-10], as well as their precision, has allowed the 
determination of PDFs across a broad range of phase space spanned by the fractional 
proton momentum carried by the struck quark, Bjorken x, and the negative squared 
four-momentum transfer between the lepton and nucleon, Q^. The high-statistics HERA 
neutral current e'^p data determine the low-x sea and gluon distributions, whereas the 
hxed-target data, taken at lower centre-of-mass energy, determine the valence distribu¬ 
tions and the higher-x sea distributions. 

The gluon PDF contributes only indirectly to the inclusive DIS cross sections. However 
it makes a direct contribution to jet cross sections through boson-gluon and quark-gluon 
scattering. Tevatron high-Fr jet data [11,12] have been used to constrain the gluon in the 
hts of MRST [1,2] and CTEQ [3]. However, these data suffer from very large correlated 
systematic uncertainties from a variety of sources. For example, the total systematic 
uncertainty of CDF data is ~ 60% over its full Et range. In the present paper, ZEUS 
neutral current e+p DIS inclusive jet cross sections [13] and direct photoproduction dijet 
cross sections [14] have been used to constrain the gluon. These cross sections have only 
~ 5% total systematic uncertainty, mainly due to the absolute energy-scale uncertainty 
of the jets. 

These jet data were used, together with ZEUS data on neutral and charged current (NC 
and CC) e^p and e~p DIS inclusive cross sections [15-20], as inputs to an NLO QCD 
DGLAP analysis in order to determine the PDFs. This ht is called the ZEUS-JETS ht. 

In the ZEUS-JETS ht, the lower NC inclusive cross-section data determine the low-x 
sea and gluon distributions^ and the high NC and CC inclusive cross sections determine 
the valence distributions. The use of ZEUS data alone eliminates the uncertainty from 
heavy-target corrections required in global analyses in which the vFe and pD hxed-target 

^ The HERA kinematics is such that the lower-Q^ data are also at low x. 


VH 



data, together with isospin-symmetry constraints between u and d in the proton and 
neutron, have been used for determining the valence distributions. The jet cross-section 
data constrain the mid- to high-x {x ~ 0.01 — 0.5) gluon PDF. The predictions for the 
jet cross sections are calculated to NLO in QCD and are used in the ht rigorously, rather 
than approximately as in previous hts [1-3]. The quality of the £t establishes that NLO 
QCD is able simultaneously to describe both inclusive cross sections and jet cross sections, 
thereby providing a compelling demonstration of QCD factorisation. 

The value of as{Mz) is hxed in most PDF hts; for the ZEUS-JETS ht, the value as{Mz) = 
0.118 [21] is used. A simultaneous ht for as{Mz) and the PDF parameters, called the 
ZEUS-JETS-a^ ht, has also been made. This ht accounts for the correlation between 
as{Mz) and the gluon shape. The addition of the jet production data provides enough 
constraints to give an accurate determination of as{Mz) despite this correlation. 

The PDFs are presented with full accounting for uncertainties from correlated systematic 
errors. Performing an analysis within a single experiment has considerable advantages 
in this respect since global hts have found signihcant tensions between diherent data 
sets [1]. In the present analysis, the contribution to the PDF uncertainties from corre¬ 
lated experimental uncertainties and normalisation uncertainties is significantly reduced 
in comparison to the previous ZEUS-S global ht analysis [8], which used data from many 
diherent DIS experiments. 

This paper is organised as follows. In Section the theoretical background is reviewed 
briehy and in Section 13 the method of analysis is outlined, with particular emphasis on 
the new features needed to include the jet cross sections in the ht. In Section HJ the 
ZEUS-JETS ht is compared to data and the extracted parton distributions and their 
experimental uncertaintes are presented. Model uncertainties are discussed and a com¬ 
parison is made to the Tevatron jet data. In Section |3 the analysis is extended to the 
evaluation of as{Mz) in the ZEUS-JETS-Os ht and the uncertainties on as{Mz) from 
theoretical sources are discussed. Section El gives a summary and conclusions. 


2 Theoretical Background 


The kinematics of deep inelastic lepton-nucleon scattering are described in terms of the 
variables Q^, Bjorken x and y, the fractional energy transfer between the lepton and 
hadron systems. The diherential cross sections for the NC DIS process are given in terms 
of structure functions by 


dxdQ'^ 


27ra^ 

xQ^ 


[Y+F,{x,Q^) 


y‘^ Fl{x, QQ =f U- xFs^x, QQ] , 
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where = 1 ± (1 — yY- The structure functions F 2 and xF^, are directly related to quark 
distributions, and their dependence, or scaling violation, is predicted by perturbative 
QCD. At ^ 1000 GeV^, the charged lepton-hadron cross section is dominated by photon 
exchange and the structure function F 2 . For x ^ 10“^, F 2 is sea-quark dominated and its 
dependence is driven by the gluon contribution, such that HERA data provide crucial 
information on both quark and gluon distributions. The longitudinal structure function 
Fl is only important at high y and is calculated, in perturbative QCD, from the quark 
and gluon distributions [22]. At high the structure function xF^ becomes increasingly 
important; it provides information on valence quark distributions. The CC interactions 
are sensitive to the flavour of the valence distributions at high x since their (LO) cross 
sections are given by 


GlM^ 


dxdQ"^ 


27rx(Q2 + m^)2 


X [(n + c) -F (1 -2/)^(d-Fs)] , 


d?a^^{e p) 
dxdQ"^ 


27ra;(Q2 + M^)2^ 


[(m-Fc) -F (1 -yf{d + s)\ 




where the parton distributions u, d, s, c are functions of x and Q^. Thus the e~p CC 
cross section gives information on the u valence quark at high x, whereas the e'^p CC 
cross section gives information on the d valence quark at high x. This is particularly 
important since this process is a direct probe of the d valence quark on a proton target at 
high Q"^. Determinations of the d valence distribution have previously been dominated by 
low data using isoscalar iron or deuterium targets. Such determinations are subject to 
uncertainties from higher-twist contributions, heavy-target and binding corrections and 
isospin-symmetry assumptions. 


The inclusive cross-section data depend directly on the quark distributions, but the gluon 
distribution affects these cross sections indirectly through the scaling violations. Pertur¬ 
bative QCD predicts the rate at which the quark distributions evolve with the scale Q 
through the DGLAP equation 


dqi{x,Q‘^) 

dlnQ2 


OisiQ^ 

2n 


dy 

y 


'^qj{y,QYP^ 


mj 


L 3 


+ a(v,Q^)P,., Cj 


( 1 ) 


where the ‘splitting function’ Pij{z) represents the probability of a parton (either quark 
or gluon) j emitting a quark i with momentum fraction z of that of the parent parton. 
Thus the gluon distribution can be obtained indirectly from the scaling violations of the 
quark distributions. The parameters that describe the gluon shape and the value of the 
strong coupling constant, as{Mz), are correlated through the DGLAP equations. 

The QCD processes that give rise to scaling violations in the inclusive cross sections, 
namely the QCD-Compton (QCDC) and boson-gluon-fusion (BGF) processes, are ob- 


IX 










served as events with distinct jets in the hnal state provided that the energy and mo¬ 
mentum transfer are large enough. The cross section for QCDC scattering depends on 
as{Mz) and the quark PDFs. For HERA kinematics, this process dominates the jet 
cross section at large scales, where the quark densities are well known from the inclusive 
cross-section data, so that the value of as{Mz) may be extracted without strong corre¬ 
lation to the shape of the gluon PDF. The cross section for the BGF process depends 
on as{Mz) and the gluon PDF so that measurements of jet cross sections also provide a 
direct determination of the gluon density. 


3 Analysis Method 

The present analysis was performed within the Standard Model conventional paradigm 
of leading-twist NLO QCD. The QCD predictions for the PDFs were obtained by solving 
the DGLAP evolution equations at NLO. These equations yield the PDFs at all values 
of provided they are parameterised as functions of x at some input scale Qo- The 
programme QcDNUM [23] was used to perform the evolution. 

The applicability of the leading-twist, NLO DGLAP formalism to HERA data was inves¬ 
tigated in the previous ZEUS analysis [8], and suitable data cuts were defined. All the 
present data lie above these cuts. The data sets fitted in this analysis and their kine¬ 
matic coverage are presented in Table d In total there are 577 data points from a total 
luminosity of 112 pb“^ from the HERA-I (1992-2000) running period. 

Full account has been taken of correlated experimental systematic uncertainties using the 
Offset method, described in the previous ZEUS-S PDF analysis [8]^. There are 22 inde¬ 
pendent sources of correlated systematic uncertainty and 4 independent normalisations 
for the data sets in the present analysis. The number of correlated systematic uncertain¬ 
ties for each data set, their normalisations and the correlations between the data sets are 
detailed in Table [TJ 


3.1 Inclusive cross-section data 

The inclusive cross-section data used in the hts were reduced double differential cross- 
sections in X and from: NG e+p scattering [15,19]; NG e~p scattering [17]; GG e+p 
scattering [16,20]; and GG e~p scattering [18]. 

^ Different treatments of experimental uncertainties in PDF analyses are discussed extensively else¬ 
where [24-26]. The Offset method gives conservative PDF uncertainty estimates. 
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The NLO QCD predictions for the structure functions, which enter into the expressions 
for the cross sections, were obtained by convoluting the PDFs with the QCD coefficient 
functions appropriate to the process. It is necessary to specify the scheme and scale 
choice for the calculations. The renormalisation and factorisation scales for the inclusive 
DIS processes were chosen to be Q. The DGLAP equations were solved in the MS 
scheme. For heavy-quark production, the general-mass variable flavour-number scheme 
of Thorne and Roberts (TRVFN) [27] was used in order to interpolate correctly between 
threshold behaviour and high-scale behaviour for heavy quarks, as discussed in the ZEUS- 
S analysis [8]. The values of the heavy quark masses used were rric = 1.35 GeV and 
rrib = 4.3 GeV. Variation of these values in the ranges 1.2 < rric < 1.5 GeV and 4.0 < 
rrib <4.6 GeV produced changes in the PDF parameters that are negligible in comparison 
to the experimental uncertainties. 


3.2 Jet data 

The jet data used in the hts were: DIS inclusive jet differential cross sections as a function 
of the transverse energy in the Breit frame, for different bins [13]; photoproduc¬ 
tion dijet cross sections as a function of the transverse energy of the most energetic jet, 
in the laboratory frame, for different jet-pseudorapidity ranges [14]. The systematic 
uncertainty from the absolute jet energy scale was fully correlated between these two sets 
of data. 

The cross-section predictions for photoproduced jets are sensitive to the choice of the 
input photon PDFs. The AFG photon PDF [28] has been used in the hts. In order 
to minimise sensitivity to this choice, the analysis has been restricted to use only the 
‘direct’ photoproduction cross sections. These are dehned by the cut > 0.75, where 
is a measure of the fraction of the photon’s momentum that enters into the hard 
scatter [14,29,30]. 

The programme of Frixione and Ridolh [31] was used to compute NLO QGD cross sections 
for photoproduced dijets and Disent [32] was used to compute NLO QGD cross sections 
for jet production in DIS. These programmes treat the heavy quarks in a massless scheme. 
However all the jet data are at scales sufficiently high that the TRVFN scheme and the 
zero-mass variable flavour number scheme (ZMVFN) are equivalent. The calculation of 
the NLO jet cross sections was too slow to be used iteratively in the ht. Thus, they were 
used to compute LO and NLO weights, d, which are independent of and the PDFs, 
and are obtained by integrating the corresponding partonic hard cross sections^ in bins of 

^ For the dijet photoproduction cross sections, the weights also included the convolution with the photon 
PDFs. 
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^ (the proton momentum fraction carried by the incoming parton), fip (the factorisation 
scale) and, for the case /iF 7 ^ fJ^R (the renormalisation scale). The NLO QCD cross 
sections, for each measured bin, were then obtained by folding these weights with the 
PDFs and a* according to the formula 

= E E E . (2) 

n a i,j,k 

where the three sums run over the order n in the flavour a of the incoming parton, 
and the indices {i,j,k) of the fip and fip bins, respectively. The PDF, fa, and 
were evaluated at the mean values (^), (/if) and (/if) of the variables /if and fip in 
each {i,j, k) bin. The factorisation scale was chosen as fip = Q for the DIS jets, and the 
renormalisation scale was chosen as fip = Ep (with fip = Q as a cross-check). For the 
photoproduced dijets, the standard scale choices were fip = fip = Ep/2 (where Ep is the 
summed transverse momenta of hnal-state partons). This procedure reproduces the NLO 
predictions to better than 0.5%. 

The predictions were multiplied by hadronisation corrections before they were used to £t 
the data. These were determined by using Monte Carlo (MC) programmes, which model 
parton hadronisation to estimate the ratio of the hadron- to parton-level cross sections 
for each bin. For the DIS jet data, an average of the values obtained using the Ariadne, 
Lepto and Herwig MC programmes was taken [13]. For the photoproduction data, 
an average of the values obtained from the Herwig and Pythia MC programmes was 
taken [14]. The hadronisation corrections are generally within a few percent of unity [13, 
14]. The predictions for DIS jet production were also corrected for contributions. 

3.3 Parameterisation of PDFs 

The PDFs for u valence, d valence, total sea, gluon and the difference between the d and 
u contributions to the sea, are each parameterised, at = 7 GeV^, by the form 

xf{x) = pix^'^{l — xy^{l + Pax). 

It was checked that no signihcant improvement in yf' results from the use of more complex 
polynomial forms or from variation of the value of Qq- The following constraints were 
imposed on the parameters pp 

• the normalisation parameters pi, for the d and u valence and for the gluon, were 
constrained by imposing the number sum-rules and momentum sum-rule, respectively; 

• the P 2 parameters, which constrain the low-x behaviour of the valence distributions, 
were set equal for u and d, since there is insufficient information to constrain any 
difference; 
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• there is also no information on the flavour structure of the light-quark sea in a ht to 
ZEUS data alone. Thus, the normalisation of the d — u distribution was hxed to be 
consistent with the measured violation of the Gottfried sum-rule [33,34] and its shape 
was hxed to be consistent with the Drell-Yan data [35]; 

• a suppression of the strange sea by a factor of two at Ql was imposed in accordance 
with neutrino induced dimuon data from CCFR-NuTeV [36,37]. 

The ht is not sensitive to reasonable variations of these assumptions, indicating that it is 
only possible to extract a havour-averaged sea distribution from these ZEUS data. 

The ZEUS inclusive cross-section data are statistics limited at large x, where the sea 
and the gluon distributions are small. This leads to sizeable uncertainties in the mid- to 
high-x sea and gluon shapes if a ht is made to inclusive cross-section data alone. The 
ZEUS jet data constrain the gluon distribution in this kinematic region. Two diherent 
strategies were used to constrain the sea distribution: hrstly, a simple parameterisation 
setting p 4 = 0 was used; secondly, the p 4 parameter was freed but the ps parameter 
was hxed to the value obtained in the ZEUS-S global ht [8]. In the latter case, model 
uncertainties on the high-x sea include the ehect of changing this hxed value of ps within 
the limits of its uncertainty as determined in the global ht. There is very little diherence 
in the shapes and uncertainties of the sea PDF as determined in these two strategies once 
this model uncertainty on p^ is taken into account. Distributions are presented for the 
former choice because of its simplicity. Finally, there are 11 free parameters describing 
the input PDF distributions, which are listed in Table |21 


4 Results 

The ZEUS-JETS ht and the NC and CC reduced cross-section data are shown in Figs.[TJ|21 
and El The ht and the jet cross-section data are illustrated in Figs. (Hand El A good 
description of the data is obtained over many orders of magnitude in scale. A measure of 
the goodness of ht for the Ohset method is obtained by re-evaluating the by adding 
the statistical, uncorrelated and correlated systematic uncertainties in quadrature [24]. 
The total obtained is 470 for 577 data points. The extracted PDF parameters and 
their experimental uncertainties are given in Table El 

The valence distributions for the ZEUS-JETS ht are shown in Fig. El Although the 
high-x valence distributions are not as well constrained as they are in global hts which 
include hxed-target data, they are competitive, particularly for the less well-known d 
valence distribution. Furthermore, they are free from uncertainties due to heavy-target 
corrections, higher-twist ehects and isospin-symmetry assumptions. 
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The gluon and sea distributions for the ZEUS-JETS £t are shown together in Fig. 0 
Whereas the sea distribution rises at low x for all Q^, the gluon distribution flattens for 
~ 2.5 GeV^ and becomes valence-like for lower Q^. The gluon and sea distributions 
are as well determined as the corresponding distributions of the global hts [1-3,8] at low 
X since the HERA inclusive NC data determine these distributions for all the hts. At 
high X, the uncertainties of the sea are constrained to be similar to those of the ZEUS-S 
global ht by the choice of parameterisation, whereas the uncertainties of the gluon have 
been reduced by the addition of the ZEUS jet data. 

In Fig. IHl the uncertainty of the gluon distribution for hts with and without the jet data 
are compared. The shapes of the PDFs are not changed signihcantly by the addition 
of jet data, even though the gluon parameterisation is sufficiently hexible to allow this, 
indicating that there is no tension between the jet data and the inclusive cross-section 
data. Although the jet data constrain the gluon mainly in the range 0.01^^ A 0.4, the 
momentum sum-rule ensures that the indirect constraint of these data is still significant 
at higher x. The decrease in the uncertainty on the gluon distribution is striking; for 
example at = 7 GeV^ and x = 0.06 the uncertainty is rednced from 17% to 10%. A 
similar decrease in uncertainty by a factor of abont two is found in this mid-x range, over 
the full range. 

In Fig. 0 the valence, sea and gluon PDFs are compared for the ZEUS-JETS ht and the 
previous ZEUS-S global PDF analysis. There is good agreement between the ZEUS PDF 
extractions. The hgure also compares the MRST and GTEQ PDFs to the ZEUS-JETS 
PDFs. These PDFs are compatible with the ZEUS PDFs, considering the size of the 
uncertainties on each of the PDF sets. 


4.1 PDF Uncertainties 

The following sources of model uncertainty have been included in the PDF uncertainty 
bands: 

• the value of Ql was varied in the range 4 < Qo < 10 GeV^; 

• the forms of the input PDF parameterisations were changed, by modihying the form 
(1 -|- p^x) to (1 -|- P 4 X -|- Psa/x) for the valence parameterisations and by varying the 
choice of constraints applied to the sea parameterisation as explained in Section El 

• the standard Et cnts applied to the jet data were raised to > 10 GeV and 
E}^^^ > 17 GeV for DIS jets and photoprodnced jets, respectively, since there are 
some small discrepancies between the ht predictions and the jet data at the lowest 
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transverse energies'^; 

• the hadronisation corrections applied to the jet data have been varied by half the dif¬ 
ference between the values obtained from the Herwig and Pythia MC programmes 
for the photoproduced jet cross sections [14] and by the variance of the values ob¬ 
tained from the Ariadne, Lepto and Herwig MC programmes for the DIS jet 
cross sections [13]. The uncertainties on the hadronisation corrections determined by 
these procedures are < 1%; they lead to uncertainties in the PDFs which are small in 
comparison to the experimental uncertainties; 

• as explained in Section HO the photoproduction data used in the £t are enriched with 
direct photon processes by the cut > 0.75; however it is not possible to select 
jet cross sections that are completely independent of photon structure. Therefore the 
sensitivity of the £t results to the input photon PDFs was investigated. In Fig. fTTIk 
the proton PDFs extracted from the ZEUS-JETS fit using the AFG photon PDFs [28] 
are compared with those extracted using the GRV [38,39] and CJK [40] photon PDFs. 
There is no visible difference in the extracted proton PDFs. In Fig.fTUb this comparison 
is shown for a fit in which the ‘resolved’ photon cross sections, < 0.75 [14], 
have been included. A significant difference is now observed between the extracted 
proton PDFs using the AFG, GRV, or GJK photon PDFs. Note that this difference 
is greatest in the region of x where the jet data have the most signihcant impact 
in reducing the uncertainty of the gluon PDF. Thus, although the addition of the 
resolved photoproduction cross sections reduces the experimental uncertainty on the 
extracted gluon PDF, it introduces a model uncertainty due to the limited knowledge 
of the photon PDFs which outweighs this advantage. Hence the present analysis used 
only the photoproduction cross sections with > 0.75. The difference in the proton 
PDFs extracted using the AFG and GRV photon PDFs was used to estimate the small 
residual model uncertainty due to the photon PDF in the ZEUS-JETS £t. 

The effect of some of the larger model variations listed above on the shapes of the extracted 
PDFs is illustrated in Fig. El These model variations are a much smaller source of 
uncertainty than the experimental uncertainties. 

In addition to these model uncertainties a variety of cross-checks have been made: 

• the minimum x of data entering the £t was raised to a; > 5 x 10“^, and the minimum 

of data entering the fit was raised to > 4.5 GeV^. These variations did not 
produce any significant changes in the PDF parameters; 

• the ZMVFN heavy quark production scheme wasused instead of the TRVFN scheme. 
The jet data are all at sufficiently high scale that the TRVFN and ZMVFN schemes 

^ This is also the case for the MRST and CTEQ PDFs [13,14]. 
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are equivalent. However, it is well known that the use of the ZMVFN scheme makes 
small differences to the shape of the gluon at a; < 10“^. This shift is well within the 
experimental uncertainty bands; 

• the choices of factorisation and renormalisation scale have been varied. The choice 

of Q is not in dispute for inclusive DIS processes. However, the scale choices for jet- 
production are not so unambiguous. Thus, factorisation and renormalisation scales 
were varied by a factor of ^2 for both the DIS jets and the photoproduced jets and 
additionally the conventional scale fir = = Q for the inclusive cross-section data 

was varied by the same factor®. The renormalisation scale for the DIS jet data was 
also changed from fiR = to fiR = Q. These changes in the choice of scale produced 
changes in the shapes of the PDFs which are small in comparison to the experimental 
uncertainties; 

• a Hessian £t was performed to the same data sets as for the ZEUS-JETS £t. The 
central valnes of the PDF parameters were fonnd to be similar to those of the ZEUS- 
JETS fit, well within the latter’s uncertainties. In the Hessian htting method [25,26], 
the theoretical prediction is used to determine the optimal correlated systematic shifts 
of the data. The correlated systematic uncertainties are assumed to be Gaussian 
distributed. This assumption is not correct for the data sets considered here, and the 
resulting uncertainties of the £t are underestimated. On the other hand, the method 
has a which is a well defined measnre of the goodness-of-ht, not available in the 
Offset method. The psr degree of freedom of the Hessian fit was 1.12 for 566 degrees 
of freedom®. 

Figure 1121 shows the ZEUS-JETS PDFs compared to those of the HI 2000 PDF anal¬ 
ysis [10]. The comparison is done in terms of the xU = x{u -1- c), xU = x{u + c), 
xD = x{d + s), xD = x{d + s) and gluon PDFs, which have been directly extracted by 
HI. The PDFs extracted by ZEUS and HI are broadly compatible. Note that the Hessian 
method of treatment of the correlated systematic nncertainties used in the HI fit results 
in a smaller experimental uncertainty on the gluon PDF [25], bnt the model uncertainty 
is significant. By contrast, the Offset method of treatment of correlated systematic un¬ 
certainties used in the ZEUS fit results in a larger experimental uncertainty, so that it 
dominates in comparison to the model uncertainties. 


® Larger variations, by a factor of 2, are not presented since they produce fits with unacceptably large 
X^- The acceptability of a is judged by the hypothesis testing criterion [8] such that the variation 
from the minimum should not exceed ~ •\/2iV, where N is the number of degrees of freedom. In the 
ZEUS-JETS fit, VM = 33. 

® If the Et cuts applied to the jet data are raised, as described in Section EH the per degree of 
freedom of the Hessian fit becomes l.OI for 554 degrees of freedom. 
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4.2 Comparison to Tevatron jet data 


It has been suggested that PDF hts to DIS data alone cannot produce a hard enough high- 
X gluon to describe the high-F' 7 ’ inclusive jet cross sections measured at the Tevatron [1], 
To investigate this issue, the ZEUS-JETS PDFs were used to make predictions for the 
CDF jet cross sections. The information on the correlated systematic uncertainties of the 
CDF data is supplied in such a way that it is possible to make a £t to these data by the 
Hessian method. In such a £t, the PDF parameters are hxed but the eight systematic 
uncertainties are freed. The of fh® CDF jet data with respect to the ZEUS-JETS £t 
was calculated using this procedure and = 48.9 was obtained. This is to be compared 
to = 46.8 which was obtained by the CDF collaboration [12], using the same procedure, 
for a £t to the CTEQ4HJ PDFs, which were specially developed to ht the CDF jet data. 
Thus, the ZEUS-JETS PDFs give an acceptable description of the CDF jet data. 


5 Extraction of as 

The strong correlation between the gluon shape and the value of as{Mz), which affects hts 
to inclusive cross-section data alone, can be broken by including the jet production cross- 
section data, which are dependent on the gluon PDF and the value as{Mz) in a different 
way from the total cross section. Jet production cross sections are directly dependent on 
the gluon PDF through the BGF process, but jet production also proceeds though the 
QCDC process, which dominates the cross section at large scales. This process depends 
on as{Mz) and the quark densities, which are directly determined from the inclusive 
cross-section data. Thus the addition of jet data allows an extraction of as{Mz) that is 
not strongly correlated to the shape of the gluon PDF. 

In previous determinations of as{Mz) using ZEUS jet data [13,41-45], the uncertainty 
from the correlation to the PDFs was taken into account by using PDFs from the global 
hts of CTEQ and MRST, which were determined assuming diherent values of as{Mz)- 
In the present analysis this correlation is directly included by htting the PDF parameters 
and as{Mz) simultaneously. The conditions for the ZEUS-JETS-cTs ht are otherwise the 
same as for the ZEUS-JETS ht. The value 

as{Mz) = 0.1183 ± 0.0007(uncorr.) ± 0.0022(corr.) ± 0.0016(norm.) ± 0.0008(model) 

was obtained, where the four uncertainties arise from the following sources: statisti¬ 
cal and other uncorrelated sources; experimental correlated systematic sources excluding 
normalisation uncertainties; normalisation uncertainties; and model uncertainty. Here the 
uncertainty on as{Mz), which usually comes from the correlation to the PDF shapes, is 
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automatically included in the experimental uncertainties. The sources of model uncer¬ 
tainty were discussed in Section IT!T1 In addition to the model uncertainties included in the 
PDF extraction, the following extra sources have been included in the model uncertainty 
for as{Mz). variation of the and x cuts on the data, as specihed in Section im and 
the use of the ZMVFN instead of the RTVFN scheme for heavy quark production. 

This extraction is at NLO. A crude estimate of the effect of terms beyond NLO can be 
made by variation of the choice of fiR. This scale was varied by a factor of y/2 for all 
the data sets entering into the £t, as described in Section im The most significant effect 
comes from the variation of the renormalisation scale for the photoproduction process. 
These scale changes produced shifts of Aas{Mz) ~ ±0.005. 

Figure El illustrates that the improved accuracy of the extraction of as{Mz) in the ZEUS- 
JETS-Q!s £t is due to the inclusion of the jet data. The prohle around the minimum 
is shown as a function of as{Mz) for the ZEUS-JETS-as £t and a similar £t in which 
the jet data are not included. The value of as{Mz) extracted is in agreement with recent 
determinations using measnrements in DIS [8, 9,13, 41-44, 46] and photoproduction of 
jets [45] at HERA and with the current world average of 0.1182 ± 0.0027 [47,48]. 

The extracted value of as{Mz) is close to the hxed value used in the ZEUS-JETS fit, and 
there are therefore no signihcant changes in the central values of the PDF parameters. 
The uncertainties of the valence and sea PDFs are also unaffected. However, there is some 
increase in the overall uncertainty of the gluon PDF because a weak correlation remains 
between as{Mz) and the gluon PDF parameters. This is illustrated for various values 
in Fig. HU 

The input of the jet data results in a much reduced uncertainty on the extracted value of 
as{Mz) compared to the previous ZEUS-a* analysis [8]. Since the present analysis was 
performed within a single experiment, the contributions from normalisation uncertainties 
and from correlated systematic uncertainties are both significantly reduced. In conse¬ 
quence, the total uncertainty on the gluon, including the uncertainty due to as{Mz), is 
reduced in comparison to the total gluon uncertainty determined in the ZEUS-Os global 
£t. 

6 Summary 

Due to the precision and kinematic coverage of the ZEUS data, it is now possible to 
extract proton PDFs and as{Mz) in a fit to data from a single experiment with minimal 
external inpnt. The ZEUS high-Q^ cross sections were used to constrain the valence PDFs, 
ZEUS low-Q^ NC data were used to constrain the low-x sea and gluon distributions and 
ZEUS data on jet production were used to constrain the mid- to high-x gluon. This 
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provides a compelling demonstration of QCD factorisation, showing that NLO QCD in 
the framework of the Standard Model is able to simultaneously describe inclusive cross 
sections and jet cross sections. The additional constraint on the gluon PDF from the jet 
production data allows an accurate extraction of the value of as{Mz) in NLO QCD, 

a,{Mz) = 0.1183 ± 0.0028(exp.) ± 0.0008(model). 

The uncertainty in as{Mz) due to terms beyond NLO has been estimated as Aas{Mz) ~ 
±0.005, by variation of the choice of scales. This is the hrst extraction of as{Mz) from 
HERA data alone. 

The total uncertainty on the gluon PDF is reduced in comparison to the ZEUS-Os global 
ht because of the greater precision of the as{Mz) measurement. The uncertainties on 
the valence PDFs are becoming competitive with those of the global hts, and they are 
not subject to uncertainties from heavy-target corrections, higher-twist contributions or 
isospin-symmetry assumptions. The precision of PDFs extracted from the global hts is 
now limited by the systematic uncertainties of the contributing experiments, whereas the 
precision of the present ht using ZEUS data only is limited by the statistical uncertainties 
and so further improvement can be expected when higher precision HERA-II data become 
available. 
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Data Set 

Ndata 

Norm 

Nsys 

Kinematic range 

of the data 

NC e+p 96-97 [15] 

242 

2% 

(1%) 

10 

1,2,3,4,5,6,7,8,9,10 

2.7 < < 30,000 GeV2 

6.3 X 10-5 <x < 0.65 

CC e+p 94-97 [16] 

29 

2% 

3 

5,6,11 

280 < < 17,000 GeV2 

0.015 < a; < 0.42 

NC e-p 98-99 [17] 

92 

1.8% 

6 

12,13,14,15,16,11 

200 < < 30,000 GeV^ 

0.005 <x< 0.65 

CC e-p 98-99 [18] 

26 

1.8% 

3 

17,18,11 

280 < < 17,000 GeV^ 

0.015 < a; < 0.42 

NC e+p 99-00 [19] 

90 

2% 

8 

12,13,14,15,19,11,20,21 

200 < < 30,000 GeV^ 

0.005 < a; < 0.65 

CC e+p 99-00 [20] 

30 

2% 

3 

17,-18,11 

280 < < 17,000 GeV^ 

0.008 <x < 0.42 

DIS jets e+p 96-97 [13] 

30 

2% 

1 

22 

125 < < 30,000 GeV^ 

8<E^ < 100 GeV 

7 P dijets 96-97 [14] 

> 0.75 

38 

2% 

1 

22 

14 < < 75 GeV 


Table 1: The number of data points (Ndata), normalisation uncertainties (Norm) 
and number of point-to-point correlated systematic uncertainties (Nsys) are detailed 
for each of the data sets used in the ZEUS-JETS fit. The kinematic regions of 
the data sets are also given. The number of independent correlated systematic 
uncertainties is specified as follows. Each independent source of uncertainty is 
assigned a number in the order of the systematic uncertainties as given in the 
corresponding publication. These numbers are given in the column headed Nsys, for 
each data set. Eor example, for the CC e'^p 94-97 data set, the first two systematic 
uncertainties are fully correlated to the fifth and sixth systematic uncertainties for 
the NC e'^p 96-97 data set. Note also that the second systematic uncertainty for the 
CC e'^p 99-00 data set is fully anti-correlated to the second systematic uncertainty 
for the CC e~p 98-99 data. The normalisation uncertainties are applied as follows. 
There are two normalisation uncertainties for the NC e'^p 96-97 data: an overall 
uncertainty and the relative uncertainty (indicated in parentheses) of the data with 
< ?>QCeV^, with respect to the higher data. The CC e'^p 94-97 data are 
dominated by the 96-97 data, so that the same overall normalisation uncertainty 
is applied to this data set. The two jet production data sets also share the overall 
normalisation uncertainty of the 96/97 data. The NC and CC e~p 98-99 data share 
a common normalisation uncertainty as do the NC and CC e'^p 99-00 data. 
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PDF 

Pi 

P2 

P3 

P4 

XUy 

(3.1±0.7±1.2) 

0.64 ±0.05 ±0.08 

4.06 ±0.18 ±0.24 

2.3 ± 1.1 ± 1.0 

Xdy 

(1.7 ±0.3 ±0.5 

0.63 ±0.05 ±0.08 

4.8 ± 0.7 ± 1.0 

2.6 ±2.2 ±2.3 

xS 

0.72 ±0.03 ±0.10 

-0.217 ±0.005 ±0.020 

7.0 ±0.8 ±2.0 

0 

xg 

(0.9 ±0.1 ±0.3) 

-0.28 ±0.02 ±0.04 

10.2 ±0.7 ±2.1 

16±4± 10 


Table 2: Table of PDF parameters at Ql = 7 GeV^, as determined from the 
ZEUS-JETS fit. The first uncertainty given originates from statistical and other 
uncorrelated sources and the second uncertainty is the additional contribution from 
correlated systematic uncertainties. The numbers in parentheses were derived from 
the fitted parameters via the number and momentum sum-rules. 
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Figure 1: ZEUS-JETS fit compared to ZEUS low-Q"^ e'^p NC reduced 
sections, ^nc- 
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Figure 2: ZEUS-JETS fit compared to ZEUS high-Q'^ NC e+p and e p reduced 
cross sections, ctnc- 
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Figure 3: ZEUS-JETS fit compared to ZEUS high-Q'^ CC e'^p and e p reduced 
cross sections, dec- 
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Figure 4: ZEUS-JETS fit compared to ZEUS DIS jet data. Each cross section 
has been multiplied by the scale factor in brackets to aid visibility. 
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Figure 5: ZEUS-JETS fit compared to photoproduced dijet data. Each cross 
section has been multiplied by the scale factor in brackets to aid visibility. 
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Figure 6: Valence PDFs extracted from the ZEUS-JETS fit. The inner cross- 
hatched error bands show the statistical and uncorrelated systematic uncertainty, 
the grey error bands show the total uncertainty including experimental correlated 
systematic uncertainties, normalisations and model uncertainty. 


XXIX 









8 

6 

4 

2 

0 

20 


10 

0 

30 

20 

10 

0 

lO '* 10^ 10"^ 10 ‘ 1 10 '* 10^ 10"^ 10 ^ 1 

X 

Figure 7: Gluon and sea PDFs extracted from the ZEUS-JETS fit. The uncor¬ 
related and total error bands are as in the caption to Eig. O 
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Figure 8: The total experimental uncertainty on the gluon PDF for the ZEUS- 
JETS fit (central error bands) compared to the total experimental uncertainty on 
the gluon PDE for a fit not incuding the jet data (outer error bands). The uncer¬ 
tainties are shown as fractional differences from the central values of the fits, for 
various values of . The total experimental uncertainty includes the statistical, 
uncorrelated and correlated systematic uncertainties and normalisations, for both 
fits. 
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Figure 9: (a) PDFs extracted from the ZEUS-JETS fit. (b) PDEs extracted 

from the ZEUS-JETS fit compared to ZEUS-S PDEs. (c) PDEs extracted from 
the ZEUS-JETS fit compared to MRST2001 PDEs. (d) PDEs extracted from the 
ZEUS-JETS fit compared to CTEQ6.1 PDEs. The total experimental uncertainty 
bands are shown for each PDE set. 
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Figure 10: (a) PDFs extracted from the ZEUS-JETS fit using different photon 

PDEs. The AEG photon PDE is used to obtain the central line, the GRV photon 
PDE gives the dashed line and the GJK photon PDE gives the dotted line, (b) 
PDEs extracted from a fit in which the resolved photoproduction cross-sections are 
included in addition to all the standard data sets for the ZEUS-JETS fit. The AEG 
photon PDE is used to obtain the central line, the GRV photon PDE gives the 
dashed line and the GJK photon PDE gives the dotted line. The total experimental 
error bands shown in these figures were obtained using the AEG photon PDE; for 
details see the caption to Eig. 0 
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Figure 11: Model variations discussed in the text are illustrated as fractional 
differences from the ZEUS-JETS central value for all the PDEs. Eor comparison, 
the shaded band shows the total experimental uncertainty. 
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Figure 12: Comparison of the PDFs extracted from the ZEUS-JETS fit with those 
extracted in the HI 2000 PDE analysis. For each analysis the total experimental 
error bands and the model error bands are included. However, the model uncertainty 
is not visible for the ZEUS-JETS fit. 
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Figure 13: The profile as a function ofas{Mz) for the ZEUS-JETS-ag fit (black 
dots) and for a similar fit not including the jet data (clear dots). The ordinate is 
given in terms of the difference between the total o^nd the minimum ~x^, for each 
fit. 
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Figure 14: Gluon distributions extracted from the ZEUS-JETS-ag fit. The un¬ 
certainties on these distributions are shown beneath each distribution as fractional 
differences from the central value. The inner error bands show the total uncertainty 
including statistical, uncorrelated and correlated experimental systematic uncertain¬ 
ties, normalisations and model uncertainties and the outer error bands show the 
additional uncertainty in the gluon coming from the variation of as{Mz). 
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